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SUMMARY 

A method recently developed for determining the steady flow of a 
nonviscous compressible fluid along a relative stream surface extending 
from hub to casing between two adjacent blades in a turbomachine is 
applied to investigate the through flow of air in an experimental 
mixed-flow impeller of high solidity. The shape of the stream surface 
is taken to be the same as that of the mean camber surface of the blade 
which consists of all radial elements. The principal equation govern- 
ing the through flow is solved by the relaxation method with the use of 
fourth-degree differentiation formulas for unequally spaced grid points 
caused by the varying hub and casing wall radii. 

A detailed analysis is made of both incompressible and compressible 
flow through the impeller, and contour plots of the stream function, 
velocity components, total enthalpy, static pressure, and Mach number 
are presented and discussed. The trends of flow variations in the 
impeller for the incompressible and compressible solutions are quite 
similar. The trend in the variation of static pressure along the casing 
obtained in the compressible solution compares very well with the experi- 
mental data obtained under the same operating condition. 


INTRODUCTION 

The axial -discharge mixed-flow compressor possesses the theoreti- 
cal potentiality of an optimum combination of high pressure ratio, 
high specific mass flow, flat operating characteristics, light weight, 
compactness, and durability. In conjunction with the experimental 
investigation of this type of compressor undertaken at the NACA Lewis 
laboratory, a method capable of analyzing accurately the nonviscous 
compressible through flow (reference l) was developed and applied to 
the theoretical flow in the experimental impeller shown in figure 1. 

A brief account of a compressible solution obtained by this method is 
given in reference 2. 
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A brief description of the method of analysis is presented herein 
considering the mean stream surface approach advanced in references 2 
and 3 , and both the inconpressible and compressible solutions are given 
in detail. The usefulness of the incompressible solution (which is 
obtainable much more quickly than the compressible solution) is evalua- 
ted and the available experimental data of pressure variation along the 
casing are compared with the theoretical solution obtained without and 
with an approximate correction for the effect of blade thickness. (For 
a more accurate correction for the effect of blade thickness on the 
through flow and for a complete three-dimensional flow analysis, the 
method given in reference 3 can be used.) 


SYMBOLS 


The following symbols are used in this report: 
a velocity of sound 

B blade thickness factor 


“B* 
n J 


gl( z ) 

H 

h 

I 

J,K,N 

M 

n 

P 

P 


differentiation coefficients in equation (12) used to multiply 
function value at point xj to give m^* 1 derivative at x-j_ 

based on n^h degree polynomial 

dcp 

dz 

1 o 

total enthalpy per unit mass, h + g- V 
static enthalpy per unit mass 
h + 2 W 2 - \ oi 2 r 2 = H - a>(V u r) 

function of r and z (equation (ll)) 

Mach number 

unit vector normal to mean stream surface 
blade pitch or spacing 
static pressure 


r,cp,z 


cylindrical coordinates relative to impeller 
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S mean stream surface 

s entropy per unit mass of air 

T temperature 

t time 

U velocity of blade , oor 

V absolute velocity of air 

W velocity of air relative to impeller 

X ratio of specific heats of air, 1.4 

p density of gas 

T blade thickness in circumferential direction 

\jr stream function defined on stream surface 

< o angular velocity of impeller 

Subscripts : 
c casing 

e exit 

h hub 

i inlet 

r,u,z radial, circumferential, and axial components, respectively 

T total, or stagnation, state 

t tip 

Superscript : 

* 


dimens ionle s s quantity 
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METHOD OF ANALYSIS 

* 

Equations Governing Flow on Mean Stream Surface 

A general idea of the average air flow through the mixed-flow 
impeller is best obtained by determining the flow from hub to casing 
along the mean stream surface, which circumferentially divides the mass 
flow in the channel formed by two adjacent blades into two equal por- 
tions (fig. 1). Because the impeller has very high solidity blades, 
for an approximate solution this mean stream surface may be assumed to 
have the same shape as that of the mean camber surface of the blade. 
With the shape of the surface known, the angular coordinate of the sur- 
face cp is then, in general, a given function of the two other cylin- 
drical coordinates r and z. However, as usual in mixed-flow and 
radial -flow impellers, the mean camber surface of the blade consists of 
all radial elements, a fact which renders cp a function of z only. 


cp = f(z) (!) 

and 

tf=Sl( z ) (2) 

(using the same notation as in reference 3), then 


W u n z 


and 


( 3 ) 


w u n u + W z n z = 0 ( 4 ) 

These four equations give the special properties for this type of flow 
surface. The turning function together with the hub and casing 

shape completely defines the impeller. 

Because the elements of the flow surface are all radial, the 
use of cylindrical coordinates results in considerable simplification 
in the determination of the flow along the surface. With the "blade 
force" in the radial direction equal to zero, the two other components 
need not be evaluated during the calculation, and the solution of the 
problem is obtained by solving the following differential equation in a 

stream function 'if, which is defined on the stream surface (see refer- 
ence 3): v 
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1 ♦ (re,)*] U - 

— 1 9r 


2~1 3 2 ^ 1 - ( r Sl) 2 8i|r a 2 -^ 


r 3r + 3z 2 


1 + (rg 1 ) : 


9 In Bp 9\|r 
&r 9r 


2J*BP It) + 2 % + <I|P £ 

az azf &J - K at 

J ar 


“ + T |i|.o 

3r ar 


(5) 


in conjunction with the relation between density and \|r- derivatives 

l 


PT,i 


I + 




1 + (rg x ) 2 


( 31 
J V 9r 




az 


r- 1 


s T,i 


-s 


% 


2 (rBp) 2 H ± 


( 6 ) 


In the preceding equations, r and z are the two independent 
variables and the bold derivative sign refers to the differentiation 
following the stream surface . The stream function \| r on the stream 
surface is related to the meridional velocity components of the air 
flowing along the stream surface by the following equations: 


W = - 

r rBp az 


(?) 


1 3+ 


W„ = 

z rBp ar 


( 8 ) 


In equations (5) to (8), B is actually related to the angular 
variation of the air velocities and the shape of the stream surface by 


In 


»i 



1 

V 


/ aw u 

\ n * w nz 



dt 


(9) 


in which the integration is performed along a streamline. In refer- 
ence 3, B is also shown to be proportional to the angular thickness 
of a stream filament whose mean surface is the stream surface con- 
sidered herein. In general, B can be determined by the method given 
in reference 3. At the start of the present investigation, however, 
this method was unavailable and B was taken to be 1. The solution 
based on B = 1 becomes exact when the ratio of circumferential 
thickness of the blade to the pitch approaches zero. In using this 
solution as an approximate solution for thick blades, the effect of 
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blade thickness should be properly corrected. (Results obtained in ref- 
erence 4 indicate that this factor B is closely related to the channel 
width ratio (P-t)/P. Inasmuch as the channel width ratio decreases 
toward 1 along the radius in the impeller, it is expected that the 
streamlines for the actual blade would everywhere rise above those 
obtained with B = 1.) 

In the solution of equation (5), the variations of I and s are 
determined by the values at the inlet to the blade and the fact that I 
and s remain constant along the streamline for adiabatic steady rela- 
tive flow of a nonviscous fluid. The values are considered to be radi- 
ally uniform at the inlet to the blade in the present solution, and con- 
sequently are uniform everywhere. The computation of density according 
to equation (6) is greatly simplified by the use of the general density 
table given in reference 3. For the simplified case of incompressible 
flow, equation (5) is solved using a constant value of density. 

Equation (5) is the principal equation to be solved for the region 
of flow where the function gq is known from the blade shape. 

Upstream and downstream of the impeller, g 1 is not known and therefore 
equation (5) is replaced by the following equation (references 1 to 3) : 


a 2 \|f 1 3i((_ 3 2 t|/ / 3 In Bp 3x|f + 3 In Bp 3\[r 

a r 2 ’ r 3r + a z 2 "l 3r 3r 3z 3z 


(rBp ) 2 
3\|f 

3r 

in which the radial variations of V u r are determined from its value 
far upstream of the impeller (equal to zero in the present case) and 
at the exit of the impeller, respectively, and from the condition that 
it remains constant along a streamline in the absence of blade action 
or of a circumferential pressure gradient. In the present solution, 
the radial derivative of V u r is zero upstream of the impeller and is 
not equal to zero downstream of the impeller. 


Wu a( y u r ) 

r ar 


?i + T^ 

dr ar 


= 0 


( 10 ) 


Method of Solution 

Equation (5) or (10) is first written in a common form 


, 3ft at + = ,, 

3r^ ar az^ 


(H) 
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The derivatives in equation (ll) are then approximated by an appro- 
priate finite -difference expression. In the present analysis, a grid 
size on the ueridional plane as shown in figure 2 is chosen, and 
fourth-degree polynomial representation is used. The finite-difference 
form of equation (ll) at grid point i (fig. 2) is 


E 


j=o 




( 12 ) 


where ^ and ^ denote the values of y on the surface correspon- 
ding to the grid points along constant-z and -r lines , respectively, on 
the meridional plane . In order to fit the grid points on the hub and 
casing walls at both ends of the impeller, it is necessary to use two 
grid sizes in the r direction (0.3544 in the inlet portion and 0.2650 
in the exit portion). At the time this investigation was begun, the 
tables of differentiation coefficients for unequal grid size given in 
reference 5 were unavailable, and the differentiation coefficients were 
computed from the formula given in that reference. 


In the present analysis, the system of finite -difference equa- 
tion (12) covering the entire domain is solved by the relaxation method 
of Southwell (reference 6) with the modification described in refer- 
ence 5. Because N in the equation is dependent on the unknown \|r, 
the solution must be improved through successive corrections. Suffi- 
ciently accurate results are obtained herein by making only three 
cycles of calculation. For a quicker and more accurate answer and 
especially for investigations where many similar cases are to be 
analyzed, it is desirable to employ the matrix method or another method 
on high-speed digital machines (references 5 and 7). Detailed steps of 
the matrix method of solution are given in reference 7. 


In the present calculations, 
less as follows: 


all quantities are rendered dimension- 



p t,i ; 


w* = 


w 


u t,i ' 


H* = 


H 


U. 


t,i 


RESULTS AND DISCUSSION 
Incompressible Solution 

The incompressible solution was first obtained for use as the 
starting value of the compressible solution and for evaluating the use- 
fulness of the incompressible solution for this type of compressor. 
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Prescribed conditions. - The variation of gp or dcp/dz directly 

computed from the mean camber surface of the experimental blade is 
shown in figure 3. It is smoothed out slightly and corrected at the 
exit for a slip factor of 0.95, which is used in the original design. 

The value of gp at the leading edge of the blade is -0.287. The 

incompressible solution is therefore good when the inlet Mach number is 
low and when 



1 

0.287 T t>i 


0.625 


or 


oo = 0.287 V z ± 

The results obtained in the incompressible solution are presented in 
figures 4 to 9. Overlays of the contour plots for the compressible 
solution discussed later are attached and may be used for comparison 
with the contour plots for the incompressible solution. 

Streamlines. - Figure 4 shows the streamlines for 10 equal divi- 
sions - of - The~maFs flow. The position of the streamlines is seen to be 
mainly determined by the shapes of the hub and casing. With respect 
to their geometrical relation at the inlet, they are first raised by 
the positive curvature of the hub wall and then lowered by the nega- 
tive curvature of the casing wall. (in this respect, it would be 
interesting to compare the streamlines in fig. 4 with the streamlines 
obtained for gp equal to zero everywhere, that is, an axially symmet- 
ric flow in the annulus with no tangential velocity.) 

Variation of velocity components . - In figure 5(a), contours of 
constant values of ^z,i are s VP er POsed on the streamlines. Because 

of the relative shift in the streamline mentioned in the preceding 
paragraph and the general decrease in channel area, the axial velocity 
of the air above the 0.3 stream l ine increases in the inducer section, 
decreases from z « 1.1 inches to z * 3.0 inches, and then increases 
toward the exit; whereas the axial velocity of the air below the 
0.3 streamline decreases in the inducer section and increases from 
there on. The air enters with a uniform axial velocity ratio of 0.625 
and leaves with a radial profile decreasing from 1.11 at the hub to 
0.91 at the casing. 

The variations of radial velocity shown in figure 5(b) resulted 
from a direct combination of the slopes of the streamlines and the var- 
iation in axial velocity. The maximum value of radial velocity ratio 
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is equal "to 0.66 along "the hub wall at a short distance downstream of 
the point of maximum slope. It then decreases rapidly along the hub 
toward the exit. 

The variation of relative tangential velocity in the impeller is 
shown in figure 5(c). At the leading edge of the blade, the velocity 
varies almost linearly with the radius because the elements of the 
stream surface are all radial and the axial velocity on the surface is 
nearly constant. This simple radial pattern is gradually distorted by 
the variations in axial velocity. At the trailing edge of the blade, 
the relative tangential velocity is nearly zero. 

The variation of the magnitude of resultant relative velocity is 
shown in figure 5(d). At the leading sdge, it increases with the radius 
because of the increase in its tangential component. Toward the exit 
where the tangential component is very small, its variation is similar 
to that of the axial component. Between the two ends, its variation is 
a combination of the variations of the two components. 

Between the 0 and the 0.4 streamlines, the magnitude of the resul- 
tant relative velocity decreases in the inducer section and then 
increases toward the exit. Between the 0.5 and the 1.0 streamlines, it 
fluctuates along the streamline. The variation of velocity along the 
casing and the hub is shown separately in figure 6. This local rise of 
velocity along the casing around z = 1.2 inches apparently is caused 
by the earlier rise of hub radius than casing radius. 

Variation of total enthalpy . - The greater change in tangential 
velocity at the tip (fig. 5(c)) overcomes the greater rise in radius at 
the hub «.nfl results in a positive radial gradient of total enthalpy with 
increasing magnitude from the inlet to z = 5 inches (fig. 7). How- 
ever, the extra turning at the lower radii between z = 5 inches and 
z = 6 inches causes the total enthalpy to have a negative radial gradi- 
ent leaving the impeller. This overcompensation in total enthalpy is 
due to the triangular trailing portion of the blade, which was designed 
for a different inlet flow. For this operating condition, a radially 
constant total enthalpy would be obtained downstream of the impeller if 
the blade terminates at, for example, the 1.4 constant total enthalpy 
line . 


Variation of static pressure. - The fact that 1 remains constant 
along a streamline for adiabatic steady relative flow of a nonviscous 
fluid gives, for incompressible flow. 


o> 2 r 2 


P + 

P 2 


2 
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constant along a streamline. At the inlet to the blade, because V u r 
is zero, this expression is equal to Hence 


or 


p _ W 2 oo 2 r 2 

p - H i " — + ~2~ 

p ~ H i r* 2 - W* 2 


In this form, the static pressure can be computed for any inlet 
values (within the range of applicability of incompressible flow) of 
H^, p, and the blade speed. Constant values of this generalized pres- 
sure parameter are shown in figure 8. Below the 0.4 streamline, the 
static pressure of air rises along its streamline when z < 2 inches 
and decreases when z > 2 inches. Above the 0.5 streamlines, the static 
pressure at first rises, drops to a minimum at z » 1 inch, and rises 
again to a maximum near z = 4.5 inches. The variation of pressure 
along the casing and hub is separately shown in figure 9. For a local 
rise of velocity along the casing at z «1.2 inches, there is a corres- 
ponding local drop of static pressure. This local drop of static pres- 
sure is undesirable in that it gives a high rate of increase in the main 
portion of the impeller for the same over-all rise in pressure across 
the impeller, and it may be the cause of the relatively low efficiency 
obtained in the experimental tests of this impeller. 

Downstream of the impeller , the pressure parameter has a positive 
radial gradient (from -0.1654 to 0.31 at casing). 


Compressible Solution 

Prescribed conditions . - In the compressible solution, an inlet 
Mach number of 0.45 is chosen, which gives a maximum relative Mach num- 
ber in the impeller of about 0.8 at the tip of the leading edge. The 
flow throughout the impeller is subsonic. For this inlet Mach number. 


P 

Pt /i 


— = 0.9055 


H, 


= 4.997 


U. 


t,i 
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The equivalent impeller tip speed at the exit is 987 feet per 
second. The first approximate solution of the compressible case gives 
streamlines having oscillatory motion of large magnitude near the leading 
edge of the impeller, which may be- caused by the sudden imposition of the 
0.287 value of gp at the leading edge of the impeller. The variation 

of gp is therefore smoothed and extended upstream according to the var- 
iation of turning along a mean streamline obtained in reference 4 (see 
fig. 18 of reference 4). This modification is shown in figure 3. 

As evidenced in the incompressible solution, overturning of the 
air may result at the hub of the trailing edge because the inlet flow 
used in this solution is different from that of the original design; it 
therefore seems best to remove the triangular section at the trailing 
edge and to determine the flow with a radial trailing edge ending at 
z = 5.5 inches. The resulting radial gradient in total enthalpy there 
will indicate how the extra turning should be added. The meridional 
section of the modified impeller used in this solution is shown in fig- 
ure 10. 

Streamlines. - Figure 10 also shows the streamlines obtained in the 
final compressible solution. In general, they lie above the incompres- 
sible ones, as required by the density term in the principal equation (5) 
or ( 10 ) . 

Variation of velocity components . - The contour plots of the velo- 
city components, the resultant velocity, and the velocity at the walls 
(figs. 11 and 12) are seen to be of shapes similar to the incompres- 
sible ones. Their magnitudes, however, are generally lower than the 
incompressible ones and the difference increases toward the exit. This 
lower value of the velocity is a combined result of the radially out- 
ward shift of streamline position and the increase in density in the 
impeller toward the exit. 

Variation of total enthalpy. - The contours of total enthalpy 
(fig. 13) again are similar to~hose obtained in the incompressible solu- 
tion. The values are higher than those for the incompressible ones, 
because a lower value of W u gives a larger change of the angular momen- 
tum of the air from its inlet value at any point in the impeller. 

Variation of static pressure . - For compressible flow, the ratio of 
static pressure to inlet total pressure is directly obtained by raising 
the corresponding density ratio to the 1.4 power. The shape of the con- 
tours shown in figure 14 is similar to that shown in figure 8 for the 
incompressible case. There is again a local drop of pressure along the 
casing near z = 1.1 inches (fig. 15) due to the local rise in velocity. 
This local drop in pressure can be avoided by increasing the casing 
radius around that location. 
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Variation of Mach number . - The contours of relative Mach number in 
the impeller are shown in figure 16. The curves are seen to be similar 
to those of the resultant velocity. For the inlet absolute Mach number 
of 0.45, the maximum relative Mach number was found to be equal to 0.805 
at the tip radius at the leading edge of the blade . The exit relative 
Mach number varies from 0.70 to 0.50 from hub to tip, corresponding to 
an absolute Mach number varying from 1.0 to 0.95. 


Usefulness of Incompressible Solution 

The similarity between the flow variations obtained for incompres- 
sible and compressible flow makes the quickly obtainable incompressible 
solution valuable in giving a good qualitative picture for the compres- 
sible case. 


COMPARISON OF THEORETICAL AND EXPERIMENTAL VARIATION 
OF STATIC PRESSURE ALONG CASING WALL 

During the experimental investigation of the mixed-flow compressor, 
the static pressure along the casing was obtained for the original impel- 
ler running at the speed and inlet condition used in the theoretical com- 
pressible solution (without the modification of an impeller at the exit 
as used in this calculation) . The variation of pressure obtained for 
two mass flows bridging the mass flow of the theoretical solution is shown 
in figure 17 . The theoretical variation obtained with the assumption 
that B is equal to 1 everywhere is also shown in the figure. The 
trends of variation of the two results is seen to be quite similar, and 
the theoretical value is seen to be higher than the experimental value 
for the same mass flow. This difference in magnitude is due partly to 
the approximation that B is equal to 1 everywhere. Actually, the 
value of B should increase from hub to casing in accordance with the 
practically uniform blade thickness and the directly increasing circum- 
ferential spacing. A more accurate solution of the present problem 
would require at least three blade -to -blade solutions of flow along three 
stream filaments of revolution from hub to casing (reference 3) . For the 
present, only an approximate correction of the variation in B is taken 
at the casing as follows : 

It is assumed that the values of 8i|r/3r and 9i|r/3z determined 
with B equal to 1 everywhere are the same as those obtained with an 
accurate variation of B, and that B is equal to (P-t)/p. The result 
is also shown in figure 17. 
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This approximate thickness correction reduces the static pressure 
at the casing less than it should because, as shown in figure 18, 
actually a^/ar is less for B = 1 everywhere than for a solution in 
which a correct variation in B is included, for with a correct vari- 
ation in B the streamline would shift toward the casing. This shift 
results in a lesser increase in W and hence a lesser drop in static 
pressure toward the exact value. Another part of the difference between 
the theoretical and experimental values is attributed to the neglect of 
viscous effects in the theoretical solution. 


CONCLUDING REMARKS 

The results obtained in the present analysis indicate that the gen- 
eral method developed for the solution of through flow in turbomachines 
offers a useful tool for analyzing the mean flow through the mixed-flow 
impeller, and that the general trend of compressible flow variations can 
be determined by an inconpressible solution, which is much easier to 
obtain than the compressible solution. 

The trend in the variation of static pressure along the casing 
agrees well with the experimental data. The local drop in static pres- 
sure along the casing at the end of the inducer section is attributed 
to the reduction in the normal distance of the stream filament of revo- 
lution at the casing, which is caused by the rise of hub radius in the 
inducer section. This local static-pressure drop can be avoided by 
increasing the casing radius near the location in question. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, January 20, 1952 


REFERENCES 

1. Wu, Chung -Hua: A General Through-Flow Theory of Fluid Flow With Sub- 

sonic or Supersonic Velocity in Turbomachines of Arbitrary Hub and 
Casing Shapes. NACA TN 2302, 1951. 

2. Wu, Chung -Hua: A General Theory of Three-Dimensional Flow with Sub- 

sonic and Supersonic Velocity in Turbomachines Having Arbitrary Hub 
and Casing Shapes. Part II. Paper No. 50-A-79, presented at the 
Ann. Meeting A. S.M.E. (New York), Nov. 27-Dec. 1, 1950. 

3. Wu, Chung -Hua: A General Theory of Three-Dimensional Flow in Sub- 

sonic and Supersonic Turbomachines of Axial-, Radial-, and Mixed- 
Flow Types. NACA TN 2604, 1952. 


14 


NACA TN 2749 


4. Wu, Chung-Hua, and Brown, Curtis A. : A Method of Analysis for Com- 

pressible Flow Past Arbitrary Turbomachine Blades On General Sur- 
face of Revolution. NACA TN 2407, 1951. 

5. Wu, Chung -Hua: Formula and Table of Coefficients for Numerical Dif- 

ferentiation with Function Values Given at Unequally Spaced Points 
and Application to Solution of Partial Differential Equations . 

NACA TN 2214, 1950. 

6. Southwell, R. V.: Relaxation Method in Theoretical Physics. 

Clarenden Press (Oxford), 1946. 

7. Wu, Chung -Hua: Matrix and Relaxation Solutions that Determine Subsonic 

Through Flow in an Axial -Flow Gas Turbine. NACA TN 2750, 1952. 


NACA TN 2749 


15 



CQ 


Figure 1. - Mean stream surface 
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Figure 2. - Meridional section of impeller and grid coordinates on meridional plane. 
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Axial distance, z 


Figure 4. - Meridional projection of streamlines obtained in incompressible solution. 
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(a) Axial component. 


Figure 5. - Meridional projections of constant velocity contours superposed on streamlines for 

incompressible solution. 
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Figure 5. - Continued. Meridional projection of constant velocity contours superposed on 

streamlines for incompressible solution. 
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(c) Tangential component. 
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Figure 5. - Continued- 


Meridional projection of constant velocity contours superposed on 
streamlines for incompressible solution. 
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(d) Resultant relative velocity. 


Figure 5. - Concluded. Meridional projection of constant velocity contours superposed on 

streamlines for incompressible solution. 
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Resultant relative velocity, W/U t 
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Figure 6. - Variation of magnitude of resultant relative velocity along casing and hub vails for 

incompressible solution. 
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Figure 7. - Meridional projection of contours of constant total enthalpy superposed on streamlines 

for incompressible solution. 
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Figure 8. - Meridional projection of contours of constant pressure function superposed on 

streamlines for incompressible solution. 
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Figure 10.- Meridional projection of streamlines obtained in compressible solution. 

of this fig. is enclosed.) 


(An overlay 
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^a; Axial componenx. 

Figure 11.- Meridional projections of constant velocity contours superposed on streamlines for 
compressible solution. (An overlay of this fig. is enclosed.) 
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(b) Radial component. 

Figure 11.- Continued. Meridional projections of constant velocity contours superposed on 
streamlines for compressible solution. (An overlay of this fig. is enclosed.) 
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(c) Tangential component. 

Figure 11.- Continued. Meridional, projections of constant velocity contours superposed on 
streamlines for compressible solution. (An overlay of this fig. is enclosed.) 
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(d) Resultant relative velocity. 


Figure 11.- Concluded. Meridional projections of constant velocity contours superposed 
streamlines for compressible solutions. (An overlay of this fig. is enclosed ' 
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Figure 12. 
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- Variation of magnitude of resultant relative velocity along casing and hub walls for 

compressible solution. 
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Figure 13.- Meridional projection of contours of constant total enthalpy superposed on streamlines 
for compressible solution. (An overlay of this fig. is enclosed.) 
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Figure 14.- Meridional projection of contours of constant pressure superposed on streamlines for 
compressible solution. (An overlay of this fig. is enclosed. ) 
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Figure 16. 


Meridional projection of contours of constant relative Mach number superposed on 
streamlines for compressible solution. 
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Figure 17. - Comparison of theoretical and experimental variation of static pressure along casing wall for compressible solution. 
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Figure 18. - Radial variation of stream function with constant 
and variable blade thickness factor. 
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